Tunneling in mesoscopic junctions 

H. O. Frota 

Departamento de Fisica-ICE, Universidade Federal do Amazonas, 
69077-000, Manaus-Am, Brazil 

February 2, 2008 



Abstract 

We have applied the Numerical Renormalization Group method to 
study a mesoscopic system consisting of two samples of metal separated 
by an insulating barrier, with nanometer dimensions, which allows the 
tunnelling of a single electron from one to the other side of the junc- 
tion. The junction is represented by a generalized orthodox model, taking 
into account the electronic scattering interaction due the hole and the 
tunnelling electrons, localized in the source and in the drain electrode, re- 
spectively. We have calculated the static properties (charge transference, 
charge average, quadractic charge average and specific heat) and the elec- 
tric conductivity of the junction, for the model parameters given by the 
tunneling matrices element Vr, the barrier energy U = /2C (where C is 
the capacitance of the system) and by the electronic scattering potentials 
Vl(r) acting on the electrons of the left(right) electrode. 



1 Introduction 



Studies of the tunnelling mesoscopic junctions in the present decade have ad- 
vanced thanks to the possibilities offered by the new manufacture technologies 
of these structures, which are increasingly lesser, heaving reached nanometric 
dimensions, resulting in capacitance lesser them 10~^^-F, which permits to study 
the effect of the quantization of the charge and electronic energy [1]. 

The orthodox model for the tunnelling mesoscopic junctions consists of two 
metal layers intercalated by a fine insulating layer. From the classic point of 
view there is no tunnelling current through the junction, due to the impediment 
of electrons to tunnel through the potential barrier represented by the insulator. 
In this condition the junction behaves as a capacitor with C capacitance . By 
applying an external potential Vext to the junction, it is loaded with a charge 
Q = CVext- This charge, which can vary continuously, is originated by the 
displacement of electrons in the electrodes in relation to the positive ions of 
metals. Even lowermost, Vext can produce small displacements of electrons, 
originating a small charge in the electrodes. The interaction between the charges 
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of these electrodes can be represented by the energy of the capacitor that is given 
by = QV2C. 

Taking in consideration the quantum effect, when the layer of insulator, in- 
tercalated enters the metal layers, will be very fine, it is possible that an electron 
tunnels from a side to the other side of the junction, through the insulator layer. 
In this case the process of electronic transport involves discrete charge of only 
one electron (e), in contrast of what occurs with the charges of the electrodes 
that are continuous. As the electrostatic energy associated with the tunnelling 
of an electron is U = e^/2C, for great values of C this energy is masked by the 
thermal energy /jbT. However, with the advance of the nanostructure manu- 
facture techniques, it is possible to construct junctions with small capacitance, 
so that the electrostatic energy U > ksT. In these conditions, the effect of the 
charge transference can be observed. Then, considering that the timnelling con- 
ductance r also is low, much smaller then the inverse of the quantum resistance 
of the junction, F <^ -Rg^, where Rq = 7rr/(2e^) ~ Q.ftKQ., the tunnelling of a 
single electron can hinder the tunnelling of the following electrons, giving rise 
to what is known as Coulomb blockade [1, 2, 3], that has excited great theoret- 
ical and experimental interest. Recent works have shown that, even at energies 
below the energy of capacitor E^, there is tunnelling due to charge fluctuation 
in the junction [4, 5, 6]. 

In the present work we have generalized the orthodox model, considering 
that the charge transference from a side (source) of the junction to the other 
(drain) creates scattering potentials in the electrodes, originated by the hole left 
in the source and by the electron that has tunnelled to the drain. The sudden 
creation of these localized potentials creates electron-hole pairs excitation at the 
Fermi level, modifying the behavior of the conductivity of the junction. The 
tunnelling process is in thermal equilibrium when the excess charge in each 
junction, originated from the tunnelling, is balanced by the phase shift of the 
conduction electrons, according to the Fridel sum rule. This problem has been 
treated similarly to the problem related to X-ray threshold singularity [7, 8], 
where the many particles states, after the timnelling, are projected on the many 
particles states before the tunnelling, what can induce X-ray infra-red divergence 
in the conductivity [8]. 

Other forms of nanostructures have deserved increasing interest from re- 
searchers, as the structure formed by two metal layers intercalated by an insu- 
lator layer, inside of which there is a metallic electrode, that usually receives 
the name of island [9]. The process of charge transport from a metal (source) 
to the other (drain) occurs by quantum tunnelling through the island, changing 
its charge by e. The efl'ect of the tunnelling of an electron through the island 
can be enhanced if the dimension of the island is very small, what confers to 
it a very reduced capacitance and, consequently, a great variation of its poten- 
tial. This high potential, that produces an electrostatic energy higher then the 
thermal energy, reduces the probability of tunnelling through the island, as it 
has been observed in granular metallic materials in the beginning of the studies 
of the tunnelling junctions [10, 11, 12], where the electronic tunnelling in low 
voltages is inhibited when the electrostatic energy of an electron of a grain is 
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greater that the thermal energy. With the development of the manufacture of 
nanostructure technology it is possible to construct metallic islands with known 
geometry, separated by tunnelling barriers, that can adequately control the tun- 
nelling of a single electron through the junction, that occurs by the jump of this 
electron from one to another island [9, 13, 14, 15]. 

Another geometry of nanostruture that has been studied is constituted by a 
metal substratum over which a layer of insulator and another one of semicon- 
ductor are deposited successively [16], with a pair of electrodes build parallel 
on the layer of the semiconductor. A positive voltage is applied to the metal 
substratum, so that the free electrons of the semiconductor arc confined to the 
superior surface of the insulator. This confinement forms a bidimensional free 
electron gas in the interface semiconductor-insulator. By applying a negative 
voltage to the pair of electrodes on the layer of the semiconductor, the move- 
ment of the electrons is limited throughout a canal in the parallel direction to 
the electrodes, thus forming a free unidimensional electron gas. Two saliences 
constructed in the electrodes have the purpose of creating potential barriers 
that confine electrons between them, producing discrete energy levels in this 
region, reason why these devices are called quantum points or artificial atoms 
[17, 18, 19, 20, 21, 22, 23, 24, 25]. 

In this work, using the Numerical Renormalization Group (NRG) formalism, 
we have developed an accurate numerical calculation for the generalized ortho- 
dox model for mcsoscopic junctions, taking into accoimt the electronic scattering 
[7, 8] . We have initiated by presenting the Hamiltonian of the Model in section 
I, where the properties that will be studied have been defined; in section II we 
have developed the formalism adopted for the diagonalization of the Hamilto- 
nian, the calculation of the static properties and the electric conductivity; the 
section III consists in the central point of this work, where the results obtained 
for the calculations of the static properties (charge transference, charge average, 
quadratic charge average and the specific heat) and of the electric conductiv- 
ity of the model are presented; in section IV we have presented the pertinent 
conchisions to this work and, finally, in the appendix we have developed a per- 
turbation theory in the tunnelling matrices element Vr for a particular case 
of the studied model, whose results are used to verify of the precision of the 
numerical calculation of the general model. 

2 The model 

The tunnelling junction studied is formed by two metals which are separated by 
an insulator layer, whose dimensions are nanometrics. The small thickness of 
the insulator allows the electrons to tunnel from a metal to the other, changing 
the charge pattern of the junction and giving origin to a tunnelling current. 
Usually the the Coulomb blockade theory assumes that the conduction elec- 
trons, before the tunnelling event, are in an equilibrium state and, after the 
timncUing, these electrons enters immediately into a new equilibrium state [1]. 
This means that the wave functions of the electronic states next to the Fermi 
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level arc not changed with the variation of the charge junction. Only a small 
displacement in the energy levels of the electrons occurs next to the Fermi level. 
This theory, therefore, fails to take in consideration the transient effect in the 
charging patterns between the two states of equilibrium. 

To account for the transient effect, some authors [7, 8, 28] have considered 
that the tunnelling event occurs in two stages, with different time scales. In the 
first stage of the electron tunnelling through the junction, the charge created 
in the surface of the metal due to the transference of an electron from one side 
to the other side of the junction, enter in equilibrium in a time of the order of 
the inverse of the plasma frequency of the metal, whose corresponding energy is 
of approximately leV. In the second stage, the electrons in the Fermi surface, 
which are associated to a response time much longer then the correspondent 
to the first stage, feel the change in the charge states of the junction as if 
this has occurred quickly and in a localized form. This difference in the time 
scale, associated with the nature of the localized potential quickly created in 
the surface of the metal, generates many particles interactions, giving origin to 
the electron-hole pairs excitations in the Fermi level. This effect is similar to 
the Mahan-Nozieres-Dominicis effect [29, 30] that occurs in the absorption and 
photoemission of X-rays in metals. In this in case an electron is ejected from a 
deep atomic level of a metal, changing the charge of an atom. The ion quickly 
created by the emission of the electron by the ray-X interacts in the final state 
with conduction electrons, resulting in a divergence in the X-ray absorption and 
photoemission spectrum. 

In this section we have presented the Hamiltonian of the model for a nanos- 
tructure mctal-insulator-metal, considering the electrostatic screen originated 
from the tunnelling of the electron through the junction. The model Hamilto- 
nian is written as 

H ^Hl + Hr + Ht + Hq + Hv (1) 

The terms Hj^ and Hr represent the left and right electrodes conduction 
Hamiltonians, respectively, given by 

HL = Y.ekclck (2) 

k 

HR = Y,e,dld, (3) 
g 

where the operators c^(cfe) and d]^{dq) create (annihilate) electrons in the re- 
spective conduction band and obey the following anticommutation relations: 
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{cfc,cfc.} = {cL4.} 



{dl,d,.} = s,,. ' ■ 

{dg,dg.} = {4,4,} = 

{Ck,dq} = 

The last of the proprieties indicates the independence of the operators of the 
two conduction bands. The term Ht describes the tunnelhng Hamiltonian 

Ht = Y^ i^Tcldq + Vfdlck) , (5) 

kq 

which allows the electrons be annihilate in the left side conduction band and be 
create in the right side conduction band, and vice versa, and Vt is the tunnelling 
matrices element, which is taken as independent of the moment for energies very 
close to de Fermi energy. The electrostatic junction energy is represented by 
the term Hq, 

H, = l^ (6) 

where C is the capacitance, Q is the junction charge, given by 

Q = l{NL-Nn) (7) 

^i=E4cfe (8) 

fe 

NR = Y,dldq, (9) 

Q 

and e is the elementary charge. 

The last term of the Hamiltonian given by the Eq.(l) corresponds to the free 
electron scattering by the charges localized in the junctions, written as 

Hv = Vl l^clck. + Vr 2^dldq., (10) 

kk' qq' 

where Vl and Vr represent the changing in the potential of the metals in the 
left and right side of the junction, respectively, when an electron tunnels trough 
the junction. 



3 The Method 



The Numerical Rcnormalization Group formalism used in this paper follows 
the method proposed originally by K. G. Wilson[31] for the solution of the 
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static proprieties of the Kondo model, and applied by Krishna-murthy, Wilkins 
and Wilson[32, 33] to calculate these proprieties in the Anderson model, later 
generalized by Prota and 01iveira[34] to obtain the dynamics proprieties of these 
models. 

The first study using the Wilson renormalization group method for nanos- 
tructures was carried through by Frota and Flensberg [5], having determined 
the fundamental state and transports proprieties for these structures, repre- 
sented by the usual model, which was not taken in consideration that when an 
electron tunnels through a junction, the conduction electrons quickly readjust 
its potential, created in the surface of the electrode, producing many particle 
effects with the creation of electron-hole pairs excitations. In the present work 
we have introduced this new ingredient and we have analyzed its effect in the 
fundamental state and transports properties. 

The method consists in the logarithmic discretization of the Hamiltonian 
of the conduction bands of the two electrodes, and in the determination of a 
new finite basis in which the Hamiltonian of the model, given by the Eq.(l), 
is written. After the iterative diagonalization, we obtain its eingenvalues and 
eigenvectors, calculate the charge average and the quadratic charge average of 
the junction, specific heat and, with the Kubo formula, wc obtain the elec- 
tric conductivity. The conduction bands are divided into logarithmic intervals 
±DK~^~^ [j = 1, 2, A > 1; 2: is a continuous parameter), and converted to a 
"hopping" Hamiltonian [32] given by 



- 2 ^ Xl^"^ (^fLf{n+l)z + (H) 



n=0 

'iJ = 2 Z^^"^ [9kz9{n+l)z + 9(n+l)z9nz J 



where Snz is obtained numerically. In this basis the electrodes are coupling 
solely via the first terms of the chains, defined as 




In the basis {,fnz,9nz}: the timncUing Hamiltonian the electrostatic 
junction energy Hq and the scattering Hamiltonian Hy are given by 
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= 2 (Vrflgo. + Vfglfo.) (15) 

/ N N \ 2 

=T E/-/--E5L5n. (16) 

\n=0 n=0 / 

(N ^ \ 

E - E 5-5- (^l/oVo. + V^iesLsoz) . (17) 

n=0 n=0 / 

Rescaling by the factor 2A~^^~^^/'^ /[{I + A^^)D] the Hamiltonians given 
by Eqs.(ll,12,15,16 and 17), the diagonahzation takes place in the reduced 
Hamiltonian Hn that is written as 

= A(^-i)/2 1 (fUin+i). + 4+1)./-) 

I n=0 

N-1 

+ E^- {s 

n=0 

+ VT{fL9oz+gUoz) (18) 



AT AT \ 2 



,n=0 n=0 
JV AT 



r, I y ^ fnzfnz ^ ] 9nz9nz 



\n=0 n=0 



where 



4W 

"-'(TOT 

The scale factor A^^^^^/^ t}jg reduced Hamiltonian is introduced to assure 
that the smallest eigenvalue is of unity order. The discrete approximation to 
the Hamiltonian of the model is given as the limit 

H= lim i(l+A-i)LlA-(^-i)/2if;v (23) 

JV— >oo 2 
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The diagonalization of the Hamiltonian Hj^i, represented by the Eq. (18), is 
carried out iteratively, using a recurrence relation which allows to calculate the 
eigenstates of the iteration (N+1) from the eigenstates of the iteration N. To 
obtain this recurrence relation, the charge operator given by Eq.(7) in the basis 
{fn-,9n} is writtten as. 



— r, I ^ ' /nz/i 
^ \n=0 




= o >] fUn. - V gLgn. (24) 



which commutes with the operator H^. since [Qn^H^] = 0. The definition of 
the charge operator is important to the numerical diagonahzation process of the 
Hamiltonian because, as Qjsi and Hiq commute, the associated matrices can 
be written as blocks of matrices that can be diagonalized independently, lowering 
the computational time. From the definition of Qat, the recurrence relation for 
the Hamiltonian iJjv is written as 

+ A^/^QN+l {nfUoz + VflffLffOz) , (25) 

where the iterative operator Jjv, which couples the iteration {N) with the iter- 
ation (A'' + 1), is given by 

In=£Nz fl,J(N+l)z + flN+l)zfNz+9^Nz9(N+l)z+9\N+l)z9Nz ■ (26) 

and g'AT is defined as 

9JV = ^ [InzInz - 9nz9nz) (27) 
From Eqs. (24 and 27), the recurrence relation for Qjv and Q% is written as 

Qn+1 = Qn + iN+i (28) 



[Qn+i? = [Qn? + 2 [Qn] Qn+i + [qN+if , (29) 

which will be useful for the calculation of the charge average and the quadratic 
charge average. 

The Eq.(25) consists in the core of the renormalization group method used 
in the present work. It permits that, once the eigenvalues and eigenvectors 
of the Hamiltonian are known, the eignvalues and the eigenvectors of the 
Hamiltonian H^+i can be determined. 

For a given iteration TV, the eigenstates consist in the following Fermi op- 
erators: /oz, fiz,---,fNz and goz, 9iz, 9Nz- Each one of these operators can 
assume two states, occupied and empty, so that for the iteration N the number 
of states is given by 2^^^+^^ and the matrices associate to the Hamiltonian Hn 
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in the basis {/„z,5«z} is of the order 2^(^+1) x 2^(^+1). For tipical values of N 
used in the present work {N = 20), the matrices that wiU be diagonahzed are of 
the order 2^^ x 2^^, which makes the numerical diagonaUzation impossible. To 
by-pass those difficulties were carried out, two procedures which arc reported as 
follows: 1) As the operator Qn commutes with the Hamiltonian Hn, then the 
eigenstates of also are eigenstates of Qn, which means that we can diago- 
nalizc in independent subspaces of Qn- Thus, the matrices associated with 
the Hamiltonian can be written in diagonal blocks that are diagonahzed 
separately, which reduces the dimensions of the matrices to be diagonahzed. 2) 
Even if diagonalizing in subspaces of same charge, the matrices continue^ with 
dimensions that makes the numerical diagonalization process inadequate. How- 
ever, as in the calculation of the static and dynamic proprieties the energies 
arc close to the Fermi energy, the procedure consists in generating the basis for 
the diagonalization of the iteration {N + 1) taking into account only the lower 
energy eigenstates of the iteration (N), according with the references [31] and 
[32]. 

The iterative process begins with the diagonalization of the Hamiltonian 
Ho, written as 



(30) 



+ 2 



\ {fLhz - 9oz9oz^ (vLfLfoz + Vngl^go, 



where only the operators foz and goz appear. In the second step, the operators 
fiz and giz are added and the Hamiltonian 



Hi = y/AHo + Io + 2UQm + Uql + qi (Vl/Uoz + Vr5oz5oz) , (31) 

is diagonahzed in relation to a basis formed by the eigenstates of the Hamiltonian 
in the first step (Hq) and the states that are constructed by the operators 
/iz' 9iz ^"^^ fiz9iz applied to the eigenstates of Hq. Following this procedure, 
the Hamiltonian H^+i is diagonahzed in relation to a basis formed by the 
eigenstates of H^ and the states given by the operator /(^jv-i-i)z' 9\n+i)z 
f{N+i)z9\N+i)z applied to the eigenstates of Hff. 

For each iteration the charge average {Qn), the quadratic charge average 
{Qn)^ the specific heat Cy, the charge transference r and the electrical current 
(t{oj) are calculated, in relation to the basis {foz,9oz,fiz,9iz, ■■■■>fNz,9Nz,}, as 
follows, 
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(Qn) 

{Q%) 

Cv 



Trcxp(-/3i?Ar) 
TtQ% cxp {-(3Hn) 



Trcxp 



N 



Tr {HNf cxp {-(3Hn) ( Ti-Hn cxp {-f^H, 



N 



Trexp 



I Trcxp 



(32) 
(33) 

(34) 



where the Boltzmann factor is written in terms of the reduced Hamiltonian 

as 



exp 



(35) 



where /3 =l/kBT', kg is the Boltzmann constant, T is the temperature and /3 a 
constant that is defined according with reference [32] , 



/3 = -/3 



-(l + A-i)A-(^-i)/2 



(36) 



Then, taking a fixed /? < 1 [32], the temperature T^, corresponding to the 
iteration A'', is given by 



T - ^ 
In — - — = 



i(l + A-i)A-(^-i)/2 . (37) 

The term Ht of the Hamiltonian H allows the tunnelling through the junc- 
tion, given by the charge transference operator r = c^dq, which is written 
in the basis {fnz,9nz} as 

T = 2,4go.. (38) 

The charge transference through the junction is given by the thermal average 
of the operator r, 



(r) = 



Tr (2,4,90,) cxp (-/?//jv) 



(39) 



Tr exp {-(3Hn) 

The electrical current operator through the junction, / , is written as a 
function of the time rate of the charge Q as 

dQ 

dt 

e d{NL - Nr) 
2 dt 



1 = 



(40) 



The above time derivative is obtained from the commutation of the operator 
Nl — Nn with the Hamiltonian H given by Eq.(l). As Nl — Nr commutes with 
Hl, Hr, Hq and Hy, then 
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--[Hi^,{Nj^-Nn)], 

and the electrical current operator is written in terms of the operators Cfe and 

dq as 

I=^Y.{^Tcidq-Vfd\ck). (41) 

k,q 

The electrical conductivity through the junction is obtained from the Kubo 
formula[29]: 

<j{'^) = -X\{F\iP)\'6{E^-E''-u:), (42) 

where uj is the frequency, and arc the fundamental and the final many 
particles states of the Hamiltonian H . with eigenvalues E^^ and E^ , respectively. 

For NRG purpose, the operator /, given by Eq.(41), is written in terms of 
the operators foz and goz (Eq. (13 and 14), respectively), 

= I (Vr^ffo. - V^glh,) , (43) 

and cr{u) (Eq.(42)) is written in terms of the eigenstates of the reduced Hamil- 
tonian Hn (Eq.(18)), as 



J2\{F{z)\I{z)\Q{z))f 
a{u;N,z) = -^5 {E^{z) - E^{z) - ^n) , (44) 

i^N / l+A-i ^-fAr-l)/2\ 

where the reduced energy Ef^{z), E^j^{z) and u)n are given by 



E^{z) = -1-^Ai^-^y^E^iz) (45) 

E^^iz) = 3-p^A(^-^)/2i;"(.) (46) 

o^N = ^^A(^-^)/^a;. (47) 

In Eq.(43) an electron is destroyed in the orbital foz of the left conduction 
band and another electron is created in the orbital qqz of the right conduction 
band, and vice versa, so that the application of the operator I{z) conserves 
the charge. In the matrices elements |(-F(z)| I{z) \^{z))\ of the Eq.(44) we take 
as initial state the fundamental state |f^(^)), and as final states \F{z)) all the 
excited states whose charges are equal to the charge of the initial state. As the 
energy levels of the conduction bands (left and right) are discretized, the lines 
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I {F{z)\ I{z) \^{z)) I that appear in Eq.(44) arc discontinuous in wjv, so that rarely 
the final energy is equal to the initial energy plus lon {E^{z) = E^{z) +lon)- If 
the numerical approach for A = 1 were possible, there would be the continuous 
limit, where the energy levels of the conduction bands would originate a dense 
spectrum for those lines. However, in this limit (A 1) the computational costs 
would be infinite. To recover the continuous limit the continuous parameter z 
[34] in the discretization process of the conduction bands was introduced. The 
continuous limit is obtained averaging in z the conductivity a{ujN,z) given by 
Eq.(44). The function S {^{z)) = 6 {Ef^{z) - E^{z) - wjv) can be written in 
terms of the roots Zi of $(^;)[36]: 



5{z - Zi) 







dz 


Z=Zi 



(48) 



The average a{ujisi) in z of the function a(ijjN,z), given by a (w at) = J^'' (j{ljOm, z)dz 
with I^z = Za — Zb, represents the continuous spectrum of the electric conduc- 
tivity, that is written as 



J2mzi)\i{zi)\n{zi))\' 

(l+AziA-(iV-l)/2j' 



dz 



(49) 



4 Results and discussion 

In this section the results of the charge transference t as a function of U, the 
charge average {Qn) as a function of Vgxt, the quadratic charge average {Q%) 
and the specific heat Cy as a function of the temperature T, and the electrical 
current a{u)) are shown. 

To verify the precision of the calculation, we present in Fig.(l) the results 

of the charge transference (r) = 2 (^fQ^Qoz^ as a function of Vt (the tunnelling 

matrices element), in the fundamental state, calculated by the NRG (black 
circles) and the exact result obtained from the Green's function equation of 
motion (full line). In the limit of Vr 0, r is null, since in this limit the 
two conduction bands are decoupled. For large Vr, the orbitals foz and goz are 
strongly coupled, disconnected from the remain of the conduction bands. In 
this case the Hamiltonian iJjv is reduced to two decoupled conduction bands 
with two energy levels, one below the base and the other above the top of the 
conduction bands, forming ligating and antiligating states, (/J^ — gl^) |0) /V2 

and (/]^ + gl^) |0) /v^, respectively, with 2 (^/d^ffoz) ^ 1- Prom Fig.(l) it can 
be noticed that the analytical and NRG results are in very good agreement, 
with an error less than 2%. 
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The results of (t) as a function of U /D calculated by NRG (black circles) 
and by second order perturbation theory in Vt (full line), for three cases: Vt = 
O.OID, Vt = O.ID and Vt = 0.3D are shown in Fig.(2). For small Vt, the 
second order perturbation theory develped in Apendices A, given by 

is in good agrcx^nicnt with NRG c;alcidation. Howc;vcr, for large Vt, the pertur- 
bation theory results are precarious, mainly for small U. 

The parameters Vl and Vr, corresponding to the scattering potentials of the 
localized charges in the mesoscopic junction, arc introduced in Fig. (3), taking 
Vt = 0.3L» and Vl = -Vr = 0, 0.1L» and 0.2D. The NRG shows that the 
scattering potentials reduce the charge transference, with more prominent effect 
for small U. 

The NRG results for (Q) as a function of eVext/U, taking Vl = Vr = 0, 
U = O.ID and Vt = O.OID are shown in Fig. (4). For zero temperature, the 
stairs, with steps spaced by U/e, is associated with the quantization of the 
charge transference in the mesoscopic junction. As Q = e{NL — Nr)/2, the 
difference between the number of particles in the two metals assumes integer 
values due the effect of the charge quantization in the tunnelling process. The 
thermal excitation smooths the stairs as is shown in the results for VbT = 0.Q2D 
(traced-point line) and k^T = 0.005£> (traced line). At high temperatures the 
thermal excitations masks the Coulomb blockade effect, reducing the spectrum 
of the charge as a function of the external potential to a straight line. The 
introduction of the scattering potential Vl = —Vr = 0.2D reduces the width of 
the steps, due the electron-hole excitations that arc created when an electron 
pass from one side to the other of the mesoscopic junction. This result suggests 
that the electron-hole excitations can reduces the Coulomb blockade, by the 
rcnormalization of the parameter U that, for these values of Vl and Vr, assumes 
the reduced value U* = 0.881i7. For Vl ^ Vr, as is shown in Fig. (6), the results 
present an asymmetry in the charge of the junction as a function of the applied 
external potential. 

In Fig. (7) is shown the quadratic charge average (Q^) as a function of T for 
Vt = O.ID, Vl = Vr = and different values of U. For fixed U, (Q^) increases 
with the temperature due the thermal fluctuation. Even for temperature zero, 
there is (Q^) originated from the quantum fluctuations in consequence of 
the tunnelling from one side to the other side of the junction. The introduction of 
the scattering potentials Vl and Vr modify the quantum fluctuation, reflecting 
in the quadratic charge average, as it is shown in Fig. (8). 

In the analysis of the charge of the junction as a function of the external 
potential we have observed that the Coulomb potential U is reduced by the effect 
of the electron- hole pairs excitation, originated from the scattering potentials Vl 
and Vr, decreasing the effect of the Coulomb blockade. From the calculation of 
the specific heat it is found out that the tunnelling matrices Vt also renormalize 
the potential U. The specific heat is shown in Figs. (9 and 10). As now we wish 
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to analysis only the effect of the tunnelling matrices Vt on the renormalization 
of U, we will consider the scattering potential Vl = Vr = 0. In Fig. (9) is shown 
the specific heat Cv/^b as a function of T for Vt = O.OID and different values of 
U, U = O.ID, 0.3-D and 0.5D. The results can be analyzed in terms of a simple 
two-level system. According to the Eqs. (6 and 7), the electrostatic energy of 
the junction is U{Nl — /A, so that when there is an excess electron, the 
energy of the junction is A = U /A. In the limit of Vt = 0, the model reduces 
to two decoupled conduction bands, separated by a energy barrier A. On the 
other hand, at the other extreme, when Vt is very large, the two conduction 
bands arc strongly coupled by the term VT{,fQz9oz + s})zf^z)^ represented 
in Eq.(18). In this limit, the orbitals /oz and c/oz are disconnected from the two 
conduction bands and the system reduces to two decoupled conduction bands, 
with two located states below the base and above the top of the conduction 
bands, identified as {Jq^ ± gQ^)/\/2, where the signals — and + represent the 
ligant and anti-ligant states, respectively. Again the system is reduced to two 
bands separated by an energy barrier A. In those two extremes, free bands and 
strongly coupled bands, the specific heat is given by the specific heat of two 
energy levels separated by A. In Fig. (9) the full line represent Cy/fcs of a two 
level system separated by an energy J7/4, which is in very good agreement with 
the NRG results for smaU Vt = 0.01£). 

In Fig. (10) Cv/ks as a function of T for a fixed U = Q.2D and different 
values of Vt = 0.2D, 0.3D, 0.4_D, 0.5D and 0.6D are shown. The results 
obtained from the NRG calculation are perfectly fitted by the specific heat of 
a two level system separated by an effective Coulomb potential U* /A. The 
renormalized U* as a function of Vt is presented in Fig. (11), which shows that 
U* ^ U vn the limit of small and high Vt, corresponding to the regimes of 
decoupled and strongly coupled conduction bands, respectively. 

The eletrical conductivity is determined by the Kubo formula and is shown 
in Fig. (12) as a function of the frequency uj, for the Coulomb energy U = O.OID, 
the scattering potential = Vr = 0, and three values of the tunnelling matrices 
element Vt = 0.05D, 0.10-D and 0.15D, with the energy scaled by U and the 
conductivity scaled by the conductivity Go for the limit of a; ^ oo [5] , given by 



The rate G/Gq tends to the unit limit as lu/U ^ 1 and drops to zero when 
lu/U 0, as is pointed out in the inset of Fig. (12). In contrast with the 
classical behavior, where there is conductivity only for w/C/ > 1, here there 
is conductivity through the junction also for energy lower then the Coulomb 
gap {(jj/U < 1), due to the quantum fluctuations that introduces uncertainty in 
the charge of the junction, as is shown in the nset of Fig. (12). We also verify 
that the conductivity is an universal function of u>/U (not shown in the figure), 
depending only on Vt- The effect of the scattering potential Vl and Vr is to 
reduce the conductivity, as it is shown in Fig. (13) for Vt = O.ID, U = O.OID 
and three different values of Vl and Vr. 



Go = 



(51) 
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5 Conclusion 



In this work we have studied the statie and dynamic properties of a mesoscopic 
junction, by using NRG to calculate the charge transference in the fundamental 
state, the charge average and the quadratic charge average, the specific heat 
and the electric conductivity of the junction. 

For Vt / D ^ 1 the orbitals fo^ and go^ are decoupled from their conduction 
bands, forming the states [fl^ — gl^) |0) /v^ and (/^^ + gl^) |0) /\/2, and the 

charge transference 2 (^fo^goz^ — * 1- For any value of Vt/D, the charge trans- 
ference is reduced by the scattering potential Vt and Vr. The charge average 
(Q) as a function of the external potential Vext presents a stairs spectrum, sug- 
gesting the charge quantization. The width of the steps is equal to U for Vl 
= = 0, and is reduced by the effect of the electron- hole pairs excitations 
for Vl = Vii ^ i). For different values of U , the results have shown that, even 
for zero temperature, quadratic charge average {Q"^} 7^ 0, as a consequence of 
the quantum fluctuations in the junction.The calculation of the specific heat 
discloses that the junction can be represented by a system of two energy levels 
separated U*. that depends on the tunnelling matrices element T. In the limit 
of Vt = or Vt oc, U* = U, since in these limits the two conduction bands 
are decoupled. For intermediate values of Vt, U* varies between zero and U. 
The electrical conductivity as a function of the frequency was obtained from 
the Kubo formula. The quantmn fluctuations originate electrical conductivity 
different of zero, even for energies lesser then U, in contrast with the classical be- 
havior, where the electrical conductivity is zero for energies lesser then U. The 
scattering potentials Vl and Vr introduce electron-hole pairs excitations near 
the Fermi level, reducing the conductivity of the junction for all frequency scales. 
The formalism developed in the present work can be extended for Hamiltoni- 
ans that represent other models of mesoscopic junction, as for instance metallic 
islands or quantum dot, with several tunnelling canals. 
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A Second order perturbation theory in Vr 



In this Appendices we have used the second order perturbation theory in Vr to 
calculate analytically the quadratic charge average and the charge transference 
in the fundamental state for the Harniltonian Hn given by Eq.(18), considering 
Vl = Vr = 0. In this case the Harniltonian Hn is written as 

Hn = H% + H^, (52) 

where 



'N-l 



I n=0 

nz9(n+l)z + 9\n+l)z9nz^ 

+ UQ% } (53) 

n=0 J 

= A(^-i)/^Ft (4,50. + glfoz) , (54) 

where is treated as a perturbation. The charge transference (r) and the 
quadratic charge average (Q^) in the fundamental state are given by 



{Q') = g^{n\H\n) (55) 
{T} = l-^{n\H\n), (56) 

where H is given by Eq.(23). 

The reduced energy of the fundamental state is 



where |r2) {\F)) is the many particles fundamental (excited) satate of -ff^, before 
(after) the tunnelling event, with reduced energ yE^j^ {E%f^). Following Wilson 
[31], the many particles states are constructed from many one body states, whose 
energy levels are given by 



% = ±A^-i for N odd (58) 

% = ±A^-^/2 for N even (59) 

In the fundamental sate |f2) all the energy levels of the two conduct bands 
below (above) the Fermi level are occupied (empty). In the final state \F), after 
the tunnelling, an electron is transfered from the level m of the left conduction 
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band, below dc Fermi level, to level i of the right conduction band, above the 
Fermi level, with energy given by 



(60) 



The operators /o,r and goz arc written in terms of the operators and 6^, which 
diagonalize the left and right conduction bands, respectively, as it follows 



k=J 



k=-J 



where 



aofe 



ao 



1- A- 



1\ 1/2 



(61) 
(62) 

(63) 
(64) 



and J = {N + l)/2. Substituting the Eqs.(61 e 62 ) into the Eq.(54), and the 
result inserted into the Eq.(57), with Epj^ given by Eq.(60), we obtain 



{n\ Hn \n) = 



'n,N 



A(^-i)/2y, 



2 J J 

EE 



=u=ir?^ + r?m + A(JV-i)/2f/ 



(65) 



To obtain H\Q) we multiply both the members of the Eq.(65) by the factor 
1±A — A~(^~^^/2, with rji and given by the Eqs.(58 and 63), respectively. 
In the limit of A — > 1, (n| is written as 



{n\H\n) = E'^-2vi 



In 



(1 + c/)^(i+^) 

From the Eqs.(55, 56 and 66) we finally have 



(66) 



-2V^ In 



-2Vt In 



Ui2 + U) 
{1 + Uf 

(2 + U) 



{2+U) 



2(l+!7) 



(67) 
(68) 
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FIGURE CAPTIONS 



Figure 1. Charge transference as a function of Vr/f, for [/ = V/, = Vr = 0. 
The full line represents the exact solution using the Green's function equation 
of motion; the black circles are the results of NRG calculation, considering the 
discretization parameter A = 3. 

Figure 2. Charge transference as a function of U/D, taking Vl = Vr = 0. 
The full lines represent the results of the second order perturbation theory in 
Vt and the black circles arc the results obtained from the NRG for A = 3. 

Figure 3. Charge transference as a function of U/D, for Vt = 0.3-D, Vl = 
— Vr = 0, Q.ID and 0.2Z), considering the discretization parameter A = 3. 

Figure 4. The stairs bahavior of the charge of the junction as a function of 
Feait, for Vt — O.OID, U — O.ID, Vl = Vr = 0, and three different temperatures, 
given in half width of the conduction banda, ksT/D. The quantization of the 
charge and the width of the steps of the stairs equal to U, which is smoothed 
by the temperature, are shown. 

Figure 5. Charge of the junction as a function of Vext, for a symetric junction 
with Vt = O.OID, U = O.ID, the scaterring potentials Vl = -Vr = 0.2D, and 
three different temperatures. The width of the steps is U* = 0.881f/, which is 
reduced by the effect of the electron-hole pairs excitation due to the scaterring 
potential. 

Figure 6. Charge transference as a function of Vext fo an assymetric junction 
with Vl = 0.2D, Vr = -O.ID, Vt = O.OID, U = O.ID and three values of the 
temperature. 

Figure 7. The square charge average (Q^) as a function of the temperature, 
for Vt = O.ID, Vl = Vr = 0, and several values of U. The values (Q^) tend 
for a finite values in the limit of T = 0, due to the quantum fluctuations in the 
fundamental state. 

Figure 8. The square charge average (Q^) as a fimction of the temperature 
T, for Vt = O.ID, U = 0.2D and several vahies of Vl and Vr. 

Figure 9. The specific heat as a function of T without Coulomb scattering, 
for small Vt and several values of U. The full line represents the specific heat 

for a two level system separated by L//4. 

Figure 10. The specific heat as a function of T for U = 0.2D and several 
values of Vt- The full line represents the specific heat for a two level system, 
separated by ?7/4 renormalized by the Vt tunnelling matrice. 

Figure 11. From the fitting of the NRG results for specific heat by the 
results of a two-level system we have verified that the capacitor energy U is 
renormalized by the tunnelling matrice Vt as is shwon in this figure. U* ^ U 
in the limits of Vr = and Vr — > oo, since the two bands are decoupled in these 
two limits. 
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Figure 12. Eletrical conductivity G scaled by Go (conductivity for large 
frequence) as a function of the frequence uj scaled by [/, without Coulomb 
scattering, with U = 0.01 and several values of Vr- 

Figure 13. Eletrical conductivity G scaled by Go (conductivity for large fre- 
quence) as a function of the frequence u scaled by U, with Coulomb scattering, 
with Vt = 0-lD, U = 0.01 and several values of Vl and Vr. 
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